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Abstract. From computational geometry comes the notion of a Gabriel graph of a point 
\m) ■ set in the plane. The Gabriel graph consists of those edges connecting two points of the 

, point set such that the circle whose diameter is the edge does not contain any point of the 

point set in its interior. We define a generalization of the Gabriel graph to n dimensions: 
the Morse poset. Using Morse theory we prove that for a generic set of 4 points in 
R 3 there are nine different Morse posets, up to combinatorial equivalence. At the end 
we mention some open questions and report on the results of computer experiments 
concerning these. We also compare our shape classification to another criterion widely 
used in computer science. 
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1. Introduction 

Take TV points Pi, • • ■ , P/v in R n and consider the function d : W 1 — > R denned by: 
in ■ d(X) = min d(X, PA 

CN) ■ 3=1,-, N 

We want to study the evolution of the sets d e = {X \ d{X) < e}, as e increases. In 
■rj- . particular we are interested in the Euler characteristic x of d e In case e is very small d e 

consists of N small solid spheres. Thus x — N. If e is very big, then d e is contractible 
J ; and hence x — 1- 

For a generic set of points, d is a topological Morse function. In that case, as e grows, d 
passes through a number of non-degenerate critical values. When d passes a critical value 
j> \ of index i, an z-cell gets attached. 

The number of critical points of index i is a^. From Morse theory we know that 

a' 



As an example, take a triangle with an obtuse angle in the plane. This is a special case 
of the above problem with n = 2 and = 3. Assume further that the two legs that 
encompass the obtuse angle have different lengths. In that case a = 3, ai = 2 and 
ct2 = 0. For an acute triangle where the edges have different lengths, we obtain a® = 3, 
ai = 3, and a 2 = 1. In this sense, there are two different generic triangles. 

Returning to the n-dimensional case, each critical point of index i corresponds to a 
subset of length i of {Pi, • • • ,P/v}, but not every subset of length i corresponds to a 
critical point of index i. ( With the obtuse triangle the 2-face of the triangle, that is: the 
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triangle itself, does not correspond to a local maximum of the function d. ) The subsets 
of length i that do correspond to critical points of d will be called active. An active subset 
defines a geometric (i — l)-simplex which we also will call active. Thus to the iV-point set 
Pi, ■ ■ ■ , P/v we can associate a set of subsets that are the active faces. This set of subsets 
is partially ordered by inclusion, and thus it is a poset. We will call it the Morse poset. 

Question: in IR n , for generic sets of N points, how many different Morse posets are 
there, up to combinatorial equivalence? 

When N = n + 1 we thus ask how many "different" generic simplices there are. This 
is a natural first problem to consider. 

In the plane the answer is two, see ||Sie99|| section 2. For obtuse triangles, the Morse 



poset is 

{{Pi},{P2h{Pa},{Pi,Pah{P2,Ps}} 
and for acute triangles, we get the Morse poset 

{{P 1 }, {P 2 }, {P 3 }, {Pi, P2}, {Pi, Ps}, {P2, Ps}, {Pi, P2, P 3 }}. 

The main theorem of this article says that in M 3 , there are nine different generic tetra- 
hedra (3-simplices). 

In the first section, we recall the relevant Morse theory. Then we establish some notation 
and state what Voronoi diagrams and Gabriel graphs are. Next, we state and prove the 
main theorem. In the last section we discuss transitions in relation to the configuration 
space of four points in IR 3 . We also report on numerical experiments concerning volume 
data of the different compartments of the configuration space where the Morse poset is of 
a certain type. Finally we compare our classification to the classification by shape types 
in IIEdeOUI. 



2. Genericity conditions 

We focus here on 4 points in IR 3 , but most of the notations and definitions have straight- 
forward extensions to the (n, N) general case. 

We write for the middle of the interval PiPj , Pijk for the center of the circumscribed 
circle of the triangle PiPjPk, and P1234 for the center of the circumscribed sphere of the 
tetrahedron. For the tetrahedron itself, that is to say the convex hull of {Pi, P 2 , P3, P4}, 
we will use the notation T. 

We impose the following 

Genericity condition 2.1. We require the set of points Pi, P2, P3, Pa to be in general 
position, so that the convex hull of P\, P2, P3, Pa is 3-dimensional. Moreover the points 
Pijk do not to lie on one of the edges of the triangle pPjPfc and also P1234 does not lie in 
one of the planes of the triangles PjPjPfe. 



The condition |2.1| means that the function d is not too badly behaved. To express more 
carefully what that means we recall the definition of a topological Morse function, see 
|Mor59|| . Let P 6 W 1 . And let / be a continuous real- valued function on W 1 . 
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Definition 2.2. / is topologically regular at P if there is some neighborhood U of P and 
a homeomorphism <fi:UQ such that one of the components of <fr is /■ The function / 
has a critical point at P G M n if / is not topologically regular at P. In that case, P is 
called a non-degenerate critical point of index i if there is a neighborhood U of P and a 
homeomorphism (f)\UO such that 



i-o /(/') ■ E 



*5 



A topological Morse function is a continuous function that has only non-degenerate critical 
points. 

For topological Morse functions the two crucial statements that hold in the differentiable 
case - the regular interval theorem and the attachment of cells, see chapter 5 in |[Mil63 



are true as well, as Morse proves in [Mor73j| . So if we show that d is topologically regular, 
we can apply those theorems, just as was done in |(Sie99 . 

We will need the following notations: 

Let 

Terr(P) = {X G M 3 |d(X,P) < d(X,P k )ior all k}, 

and Vi = Terr(Pj), Vy = V$ fl Vj (i different from j), V^k = V, D Vj fl (all three i,j,k 
different), V 12U = V x n V 2 D V 3 n V A . 



Proposition 2.3. The function d is a topological Morse function if the condition is ful- 
filled. In that case, d is topologically regular in all points of ofM. 3 , except in P\, P 2 , P3, P4, 
where d has a minimum and (perhaps) in the points Pij, Pijk and P1234 ■ Moreover d has 

- a minimum exactly in the points Pi, P2, P3, P4 

- a 1-saddle (saddle point of index 1) in Pij iff. Pij = fl PiPj 

- a 2-saddle (saddle point of index 2) in P^k iff. Pijk = Vijk D PiPjPk 

- a maximum in P1234 iff- P1234 = ^1234 H T, equivalently P1234 G T. 

Proof. If x = Pi then because the points lie in general position all Pj with i ^ j lie at 
some positive distance from x, so d has a minimum there. If x G Terr(Pj) but x 7^ Pj 
then x is obviously topologically regular. 

The function d restricted to the interior Vij has a minimum if P^ lies in that interior. 



This can only happen when P^ 7^ Pijk, which is assured by |2.1| . 

In the directions, orthogonal to V^, d decreases so we see that d has a critical point of 



index 1 at P^. 



The function d restricted to the interior of V^ has a minimum at P, L jk if Pijk lies in 
that interior. This can only be the case when P^ 7^ Pi234- In the directions orthogonal 
to V^k d decreases, so d has a critical point of index 2 at Pijk- 

Finally if P1234 = V1234, d obviously has a maximum there. □ 



Remark 2.4. Closer inspection might yield that |2.1| is actually not necessary. See ||Sie9£ 
for the 2-dimensional case. 
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Throughout this article we will always assume that point sets satisfy the genericity 
condition 127 



If d has a critical point as described by the conditions of the proposition ^3] then we 
say that the corresponding center is active. As explained in the introduction, an active 
center determines a subset of the points Pi, P2, P3, P4, which we call active subset. The 
Morse poset is the set of active subsets. Two sets of points in R 3 are called combinatorially 
equivalent if there exists a bijection, that sends the active subsets onto each other. We 
want to give a classification with respect to this equivalence relation. 

3. Morse theoretic possibilities 



From proposition |0| we know the maximal number the critical points of each type. 
Moreover the Euler characteristic should be +1. This gives a priori the following 9 pos- 
sibilities: 



m 


si 


s2 


M 


4 


6 


4 


1 


4 


5 


3 


1 


4 


4 


2 


1 


4 


3 


1 


1 


4 


2 





1 


4 


6 


3 





4 


5 


2 





4 


4 


1 





4 


3 









But not all possibilities will occur. Since we start with 4 points and the result should be a 
connected space, we need at least 3 saddle points of index 1. This rules out the possibility 
(4,2,0,1). 

We remind the reader that the definitions that follow assume the genericity condition 
|2.1| . For the more general definitions see | |BKUS97| ] or | [Ede01| 



The Voronoi tesselation of a point set {Pi, • • • , P/v} in W 1 consists of the union of the 
sets Terr(Pj), together with their natural combinatorial structure. Its (n — l)-skeleton 
consists of the intersections Terr(Pj) nTerr(Pj) and is usually called the Voronoi diagram 
The Delaunay "triangulation" of that point set is the simplicial complex dual to the 
Voronoi tesselation. Its 1-skeleton consists of the line segments 



{PiPj 

The Gabriel graph ( see 



dim(Terr(Pj) D Terr(Pj)) = n — 1 } 
BKOS97| , where it is defined in M? ) of a point set in 



is 

formed by using the points as vertices and placing an edge between two vertices exactly 
when the sphere, whose diameter is given by that edge, does not include any point from 
the set. The Gabriel graph is a subset of the Delaunay triangulation. Clearly the Gabriel 
graph is always connected. 

Our Morse poset is related to these much-used notions from computational geometry 
as follows. If in our case the edge between Pj and Pj is part of the Gabriel graph then, 
the point P^ is a saddle point of d. The Morse poset is a generalization of the Gabriel 
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graph. The Gabriel graph consists of the active subsets of length 2. Two Morse posets 
can only be combinatorially equivalent if the underlying Gabriel graphs are the same. 

4. Activity conditions 

We list in figure [l] the (a priori) possible Gabriel graphs for the above cases, they are 
the connected graphs with 6 vertices. Just as the case (4, 2, 0, 1) can not occur, we will 




Figure 1 . List of graphs 

prove in this section by deriving a contradiction that the case (4, 4, 2, 1) P, (4, 3, 1, 1) T 
and (4, 3, 1, 1) L do not occur. We will prove our main theorem: 

Theorem 4.1. Up to combinatorial equivalence of their Morse posets there are nine 
generic tetrahedra. They are uniquely described by the nine Gabriel graphs (4,6,4,1), 
(4,6,3,0), (4,5,3,1), (4,5,2,0), (4,4,2,1)0, (4,4,1,0)0, (4,4,1,0)P, (4, 3, 0, 0) L and 
(4, 3, 0, 0) T, drawn in figure [J. 

Proof. In the three cases that are to be excluded P1P2P3P4 and P1P2P3 are active. Let 
us see what these two conditions mean. 

4.1. Saddle points of index 1. We consider the plane E through P1P2P3 (and put is as 
'ground plane' in the picture). The half space containing P4 is called 'above' , the other 
is called 'below'. 

We are going to consider the condition that the point Y = P i4 is active. We first make no 
assumptions about the triangle P1P2P3, except that it does not have a right angle. We 
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consider the point X = P 4 as a variable and denote its projection on the P-plane by X'. 
Let Y = P14 and its projection Y'. From the condition that Y = P14 is active we get: 

d(P 1 , Y) = d(P 4 , Y) < min {d(P 2 , Y), d(P 3 , Y)} 

For the projection, this means that in terms of Voronoi diagrams in E, 

Y' E Terr (Pi) 

In terms of the projection X', this means: 

X' G 2 Terr (Pi), 

where with 2Terr(Pi) we mean scalar multiplication of Terr(Pi) by a factor 2 from the 
point Pi. 

We can do the same for the activity of P24 and P34. We get three subsets 2Terr(Pi), 
2 Terr(P 2 ), 2 Terr(P 3 ), which cover the plane E. They divide the plane into regions, where 
one, two or three of the points P14, P 24 or P 34 are active. 

Inside the plane E, the border of Terr(Pj) consists of two half-lines that meet in P123. 
The scalar multiplication by two maps P123 to P*, the antipodal point of P, on the circle 
through Pi, P 2 and P3. 

Next, look at the case of an acute triangle drawn in figure |[ 

The regions where only one P^ is active are both outside the disc D, which is bounded 
by the circumscribed circle of triangle P1P2P3. 
The picture for the obtuse case is in figure ^. 




Figure 2. The plane E when the triangle is acute 
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Pi 4 active 




Figure 3. The plane E when the triangle is obtuse 

Let P3 the obtuse angle. We see that the region where only P34 is active is outside 
D, but that the regions where P14 or P 24 is the only active saddle point can have some 
intersections with D. 

4.2. The center M of the circumscribed sphere. We assume again that P1P2P3 lies 
in the plane E. We know that M = P1234 belongs to the axis of the triangle P1P2P3, 
so its projection is M' = Pi 23 . Fix M for the moment, and consider X = P 4 as a 
variable. Possible positions of X = P4 are on the sphere with center M and radius 
r = rf(Pi,M) = d(P 2 ,M) = d(P 3 ,M). Denote by P^P^P^, the mirror image of PiP 2 P 3 
after central reflection through Pi23- 

In the cases we wish to exclude, M is active, i.e. M is a point, where the distance 
function has its maximum, i.e. M lies inside the tetrahedron. The point X must lie 
in the cone with top M and base the triangle P1P2P3, but also on the sphere (M, r). 
Consider the the plane E" , which is the image of E by central reflection through M. E" 
is parallel to E. The cone intersects the plane in a triangle, which projects (orthogonally) 
exactly onto P 1 *P 2 *P 3 *. X belongs to a part of the cap of the sphere, which lies "above" E 
(seen from M). It follow that X projects to a point X' which lies inside the disc D. (NB. 
the projection area contains at least the triangle P 1 *P 2 *P 3 *.) This is true for any position 
ofM. 

4.3. Certain cases don't exist. For the acute case this is all we need. In the obtuse 
case though we can be more precise. Let P3 be the obtuse angle. It is necessary for the 
activity of M that X' belongs to the circle sector P 2 *M'P 1 * (the one which does not contain 
P 3 ). We combine the activity conditions for M and one single P i4 . 
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Start with an acute triangle in the ground plane. It follows that the combination M 
active and a single P4 active can not occur together. 

This rules out the graphs: (4, 3, 1, 1) T , (4, 3, 1, 1) L and (4, 4, 2, 1) P. 

4.4. Positions of index 2 saddles. We are left with nine possibilities of Gabriel graphs. 
The Gabriel graphs do not give a complete picture of the combinatorics of the Morse 
points. The Gabriel graph, together with the information about activity of M, tells us 
only the number of saddles of index 2, but not the position. We discuss this now in detail 
and show that there are in all cases unique positions for the saddles. 
(4,6,4,1) : unique places (no choices). 

(4, 5, 3, 1) : There are two triangles, where all three midpoints are active. Both triangles 
must be acute. Take one of them in the 'ground plane' and assume it to be P1P2P3. Since 
P1234 is active, P1234 lies on the same side as P4. It follows that V123 must intersect the 
ground plane in the point P123, so that point is active. The same reasoning applies to the 
second triangle; this fixes the place of two 2-saddles. The two positions left for the third 
saddle are combinatorially equivalent. 

(4. 4. 1.0) : all places equivalent. 

(4.4.2.1) : Suppose P i2 is not active. Now P 3 or P 4 must be contained in B = 
B(Pi2,ru), where Tvi = d(Pi,Pi 2 ). If P123 is active, then triangle P1P2P3 must be acute, 
this means that P 3 is outside the ball B. If P 12 4 is active, then triangle P 1 P 2 P4 must be 
acute; this means that P4 is outside the ball B. It follows that the situation where two 
active 2-saddles are "separated" by an non-active edge cannot occur. This fixes the places 
of the 2-saddles up to permutation. 

(4, 6, 3, 0) : Since all midpoints of ribs are active, we have that all triangles are acute. 
Since P1234 is not active it lies outside T. There is exactly 1 triangle such that P1234 and 
P4 lie on different sides of the plane of that triangle. The corresponding P^ on that axis 
cannot be active. This fixes the places of the other 2-saddles. 

(4,5,2,0) : Choose a triangle, say PiP 2 P 3 , where all three edges are active. If P123 is 
not active, then it follows that P1234 lies below the ground plane E. We look again at 
the projection P' 4 . First P' A must lie inside D. But since two of the Pj4 (say P14 and P24) 
must be active we know that P' 4 also must lie in the activity region, described above. The 
intersection is a subset of D, which is contained in the region D* bounded by the arc 
PiP 3 *P 2 and the interval P2P1. There are still other conditions to meet: 

- P 4 must lie outside the ball P(Pi 2 ,ri 2 ), since Pi 2 is active. 

- P 4 must lie inside the ball P(P 12 3, T123), where r 12 3 = rf(Pi,Pi 2 3). 

This is not simultaneously possible if P' A lies in D*. It follows that there is only one 
possibility: the 2-saddles are the centers of two triangles with all edges active. 

(4, 4, 1, 0) P : Let PiP 2 P 3 be the triangle with all three midpoints of the edges active. 
This triangle is acute, take its plane as ground plane. If P1234 is above the ground plane, 
then P123 is active. Assume now P1234 is below the ground plane. It follows that P' A lies 
inside the disc D. But the fact that only P 14 is active means that P' A is outside. This is a 
contradiction, so only P 12 3 can be active. 

(4, 3, 0, 0) T : all places equivalent. 
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(4, 3, 0, 0) L : all places equivalent. 

That all these cases do occur follows from the computer experiments described below. 



The proof of theorem |4.1| is complete. □ 

5. Notes and remarks 

5.1. Higher dimensional results. We have not been able to prove a classification the- 
orem in ¥L n . However, upon request the authors will send interested readers a computer 
program that calculates the list in higher dimensions by just trying a lot of random point 
sets; see below. 

Except for the results in [ |Sic99|| on 4 points in the plane we have no results on the 



number of Morse posets for N points when N > n + 1. 

5.2. Relation with computer science criteria. In the above it was remarked that the 
Morse poset is a generalization of the Gabriel graph to higher dimensions. Our notion of 
a Morse poset, being new, has not yet an efficient stable algorithm to compute it for a 
general point set. 

5.3. The configuration space. We consider the configuration space of 4 points in M 3 . 
This is M 3 without the generalized diagonal 

Theorem 5.1. The Morse poset does not change under scaling, rotation, or translation. 
The quotient of (R 3 )^ by these group actions is a smooth 5 -dimensional space. 

The proof is not very difficult, and we refrain from stating it here. The non-generic 
tetrahedra form a hypersurface in the configuration space. The theorem above says that 
there are nine different types of compartments in the complement. It does not give us 
any information on the number of components of the complement and their topology. 

We do have the following information on the adjacency of the types of compartments, 
for general n and N = n + 1. 

Theorem 5.2. Let {P(X)i} i= i t ... n+1 be a generic smooth path in the configuration space, 
that intersects the hypersurface exactly once in the interval [0,1]. // the Morse poset 
changes there is a j such that the Morse poset of {P(l)j}j = i j ... jn+ i has one more ( resp. 
one less ) active subset of length j and one more ( resp. one less ) active subset j + 1 than 
{/'(('), 1, ,....„.,. 

Again the proof is not very difficult, so we do not state it here. 

What the theorem says for instance is that the two types (4, 3, 0, 0) with different 
Gabriel graphs are not adjacent. 

5.4. Statistics. We carried out some statistical experiments to see how the different types 
of tetrahedra are distributed among different 4-tuples of points. In this way complement 
of the discriminant volume data are obtained. Four points in R 3 do not form a bounded 
space, and thus there lives no uniform probability distribution on it. As the Morse poset 
is invariant under translations, scaling and rotations, without loss of generality we can 
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assume that all 4 points lie on S 2 . By translation and scaling, this can always be achieved. 
The ratio of the ( infinite ) volumes of the complement of the discriminant space will not 
change. 

A reference for the following material is [|Spi79|1 , chapter 9. 

In general, if one has a Riemannian manifold (M, g) the metric induces a volume form 
on M. The simplest Riemannian metric comes when M is embedded by a map 7 as an 
orientable hypersurface in M. n . Then the volume is given by: 



det (gij) dx = / ||7*ei A • • • A 7*e n _i|| 2 dx 
>M ' Jm 

Suppose we have a map (M/Z) n_1 — ► W 1 that is onto Mcl 
on [0, l] n_1 leads to a uniform distribution on M iff. ||7*ei A 
function on M. Let 



The uniform distribution 



A 7*e 



71—1 



2 is a constant 



7(01, a 2 ) = (sin(arccos(2ai — 1)) sin(27ra2), sin(arccos(2ai — 1)) cos(27ra 2 ), 2ai — 1) 

With this choice it turns out that dV = Andai A da 2 . So the induced probability measure 
on S 2 is uniform. We take two random numbers in [0, 1] and map them to S 2 using 7. 

For our experiment we used the Gnu Scientific Library, see | |G + 04| | . This library has an 
implementation of the apparently very reliable MT19937 random number generator. We 
took samples of 10 8 tetrahedra. Here is one. 



1 


(4,3,0,0) 


L 


17,807,919 


2 


(4,3,0,0) 


T 


898,689 


3 


(4,4,1,0) 


O 


26,224,574 


4 


(4,4,1,0) 


P 


16,421,773 


5 


(4,5,2,0) 




24,350,101 


7 


(4,4,2,1) 




3,266,345 


6 


(4,6,3,0) 




1,797,721 


8 


(4,5,3,1) 




2,697,783 


9 


(4,6,4,1) 




6,535,095 



Other samples gave approximately the same numbers, with a maximum difference of 3000. 

For three random points on the circle, the chances are 50 percent that one gets an 
obtuse triangle. So these results are very different, and we have no explanation for them. 



5.5. Edelsbrunner ratio. 

cumsphere of T is R. 



Denote, for brevity, dij = d(P{,Pj). The radius of the cir- 



Definition 5.3. The Edelsbrunner ratio p is the circumradius R divided by the minimal 
edge length min d^ . 



The ratio is used by Edelsbrunner ( see |[Ede01|| , section 6.2 ) to classify tetrahedra into 
"shape types". This article has the same objective, so it is worthwhile to compare his 
criterion to ours. We will show that small values of p can be attained by all the types 
listed in theorem |4.1|. Depending on one's taste one can conclude that our classification is 
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finer than the one proposed by Edelsbrunner, or that it describes other features. In any 
case, we hope to have improved on what Edelsbrunner calls a "fuzzy undertaking" . 
Up to rotations and translations, the simplex is determined by the lengths of its six 



edges. We assume that the simplices satisfy the genericity condition |2.1| . The generic 
simplices are an open subset V of (lR>o) 6 . 

Lemma 5.4. Consider the circumradius R as a function of the lengths d^. It is defined 
when the four points are not coplanar. It is homogeneous of degree 1 and, where defined, 
it has nonzero gradient. In addition, on V, we have 

QD 

(5.1) 7^r>° 1 < * < J < 4 

Proof. The first two statements are obvious. For the gradient note that 

d{j ~jrz = R{d\2i ' ' ' i d^i) 

i<i<j<4 m v 

because R is homogeneous of degree one. The circumradius is always > 0, so the gradient 
is never zero. To see ( |5.1|) , fix the position of Pi, P2 and P3. Varying only 0^34 the point 
P4 moves along a circle that has axis P1P2. Along this circle, R has a "maximum" when 
the four points are coplanar, and a "minimum" when P3 = P4. Both cases correspond to 
non-generic simplices. In between, the partial derivative of R wrt. 6^34 is clearly > 0. □ 

As a consequence, the Edelsbrunner ratio is a homogeneous function of degree 0. Hence, 
it is no restriction to assume R = 1. The above lemma says that {R = 1} is a smooth 
manifold W in V. The following lemma can also be found in [|Ede01| . 

Lemma 5.5. p has a unique minimum on W when all the dij are equal. 



Proof. Note first that it follows from ( |5.1|) that the d^ are regular functions in the sense 
of ^]2|. We thus have six regular functions on a five dimensional space of which p is the 
maximum: 

p = maxf— ... 1 ^ 1 

\dv2 'd 34 / min(di2, • • • , du) 

Suppose that, on an open set V in R N , we have N + 1 regular functions / , • • • , /at. Let 
/ be the maximum: / = max(/ , • • • , /at). We assert that the function / can only have 
a minimum at x 6 V when f(x) = fi(x) , for % = 0, • • • , N. 

Indeed because the functions are regular we can assume that they are all linear func- 
tions. Our assertion is trivial in that case. 

It follows that p can only attain its minimum when all d^ are equal. □ 

Lemma 5.6. The minimal edge length is always the length of an active edge. 

Proof. Suppose the minimal edge length were achieved by a non-active edge, with end- 
points X and Y. The edge is not active, so inside the circumsphere corresponding to that 
edge there is a point that is not an endpoint of the edge. However the distance from that 
point to either one of the endpoints X and Y is smaller than twice the circumradius. □ 
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Theorem 5.7. On each of the compartments of the configuration space corresponding 
to theorem \4-l\ p is bounded from below by the values in the table below. The infimum 
corresponds to the quadruple of points in the third column. These quadruples lie on the 
discriminant hyper surf ace from theorem \5. except for the case 4641, which corresponds 
to the global minimum of p. 



Type 



(4,3,0,0)L 
(4,3,0,0)T 

(4,4,1,0)P 

(4,4,1,0)0 
(4,5,2,0) 
(4,6,3,0) 
(4,4,2,1) 
(4,5,3,1) 
(4,6,4,1) 



For the case 



Rho 



§73 
|V3 

12 

lv/6 



Infimum 



(0,0,0), (1,0,0), (1,1,0), (1,1,1) 
(1,0,0) , (0,1,0) , (0,0,1), (0,0,0) 

(\y/?> -\, o),(-\Vz-\), (0,1,0), (0,1,V3) 

(0,0,0), (1,0,0), (1,1,0), (0,1,0) 

(0,0,0), (1,0,0), (1,1,0), (0,1,0) 

( cos a j, sin a j,0) j = 1, • • • , 3, (0,0,1) 

(1,0,0) , (0,1,0), (0,-1,0) , (cosa,0,sina) 

(1,0,0) , (0,1,0), (0,-1,0) , (0,0,1) 

(0, 1, 0), (lV3,-iO),(-lV3,-l, 0), (0, 0, y/2) 



'4,6,3,0) the triangle P1P2P3 is acute. For the case (4,4,2,1) the angle 



a should satisfy — ~ < a < 0. 

Proof. The case (4, 6, 4, 1) has been handled in the above. 

Because the minimum is unique, the other cases only have an infimum for p. These 
infima should thus correspond to quadruples on the discriminant hypersurface. 

The cases where M = are characterized by P1234 ^ CH(Pi, • ■ ■ , P4). We may assume 
that P1234 is the origin and that R = 1. It readily follows that all edges have length > y/2. 
For the cases (4, 4, 1, 0) O, (4, 5, 2, 0) and (4, 6, 3, 0) the configurations in the above table 
show that this is actually the infimum. 

If there is one triangle with an obtuse angle, we obtain in the same way that p > \y/2. 
Thus, the configurations for (4, 4, 2, 1) and (4, 5, 3, 1) in the above table are actually infima. 

We now need to study the first three cases. Let us start with (4, 4, 1,0) P. We will 
assume that the triangle P1P2P3 is active as well as the edges P1P2, P1P3 , P1P4 and 
P2P3. Maximizing the minimal length of the four edges on the Gabriel graph results in 
four edges of equal length. Take an equilateral triangle in the plane: 

P l = (0, 1, 0) P 2 = (-~V3, -i 0) P 3 = (iV3, -~ 0) 

This triangle has sides of length a/3. Drawing the figure in this case gives figure (|. Thus 
the projection of P4 to the P1P2P3 plane should lie in the gray area above P\ in figure |j. 

Hence all cases with the projection of P4 in the colored region andd(Pi,P 4 ) = give 
(4, 4, 1,0) P. The infimum occurs when the circumradius is minimal, and this is when 
P 4 = (0,1,^)- In that case 

P= {T2 

The (4, 3, 0, 0) cases. There are two cases: either none of the triangles are acute, or at 
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Figure 4. 4410P Plane 



least one of the triangles is acute. 

If none of the triangles are acute we look at figure |3[ It follows that the projection of 
P4 should lie outside of the triangle P1P2P3. Hence the active edges are P1P2, -P2-P3 and 
P3P4. To maximize the minimal edge length of these active edges, all three should have 
the same length, and the triangles P1P2P3, P2P3P4 should be right angled. The simplex 
is then fixed by the values in the table above. 

In the second case, assume P1P2P3 is acute. The triangle is not active, so P4 lies inside 
the sphere of P1P2P3. Look at figure The Gabriel graph contains the smallest edges, 
so all the other triangles are obtuse angled. So the active edges are P1P4, P2P4 and P3P4. 
To maximize the minimal length of these edges, they should all have equal length. We 
see that the infimum for p is as in the table above. □ 

Remark 5.8. With a little heuristic reasoning, omitted here, it can be argued that the 
configurations in the above table are actually unique infima. 
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